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ABSTRACT 



We have obtained millimeter wavelength photometry, high-resolution opti- 
cal spectroscopy and adaptive optics near-infrared imaging for a sample of 26 
Spitzer-selected transition circumstellar disks. All of our targets are located in 
the Ophiuchus molecular cloud (d ~125 pc) and have Spectral Energy Distri- 
butions (SEDs) suggesting the presence of inner opacity holes. We use these 
ground-based data to estimate the disk mass, multiplicity, and accretion rate for 
each object in our sample in order to investigate the mechanisms potentially re- 
sponsible for their inner holes. We find that transition disks are a heterogeneous 
group of objects, with disk masses ranging from < 0.6 to 40 Mjup and accretion 
rates ranging from <10~^^ to 10~^ Mgyr"^, but most tend to have much lower 
masses and accretion rates than "full disks" (i.e., disks without opacity holes). 
Eight of our targets have stellar companions: 6 of them are binaries and the other 
2 are triple systems. In four cases, the stellar companions are close enough to 
suspect they are responsible for the inferred inner holes. We find that 9 of our 
26 targets have low disk mass (< 2.5 Mjup) and negligible accretion (< 10~^^ 
MQyr"^), and are thus consistent with photoevaporating (or photoevaporated) 
disks. Four of these 9 non-accreting objects have fractional disk luminosities < 
10~^ and could already be in a debris disk stage. Seventeen of our transition 
disks are accreting. Thirteen of these accreting objects are consistent with grain 
growth. The remaining 4 accreting objects have SEDs suggesting the presence of 
sharp inner holes, and thus are excellent candidates for harboring giant planets. 

Subject headings: circumstellar matter — binaries: general — planetary systems: 
protoplanetary disks — stars: pre-main sequence 
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1. Introduction 

Multi-wavelength observations of nearby star-forming regions have shown that the vast 
majority of pre-main-sequence (PMS) stars are either accreting Classical T Tauri Stars 
(CTTSs) with excess emission extending all the way from the near-IR to the millimeter 
or, more evolved, non-accreting. Weak-line T Tauri Stars (WTTSs) with bare stellar 
photospheres. The fact that very few objects lacking near-IR excess show mid-IR or 
(sub) millimeter excess emission implies that, once the inner disk dissipates, the entire 
primordial disk disappears very rapidly (Wolk & Walter 1996; Andrews & Williams 
2005,2007; Cieza et al. 2007). The few objects that are caught in the short transition 
between typical CTTSs and disk-less WTTSs usually have optically thin or non-existent 
inner disks and optically thick outer disks (i.e., they have reduced opacity in the inner 
regions of the disk). 

The reduced opacity in the inner disk is the defining characteristic of the so called 
transition disks. The precise definitions of what constitutes a transition object found in 
the disk evolution literature are, however, far from homogeneous (see Evans et al. 2009 
for a detailed description of the transition disk nomenclature). Transition disks have been 
defined as objects with no detectable near-IR excess, steep slopes in the mid-IR, and large 
far-IR excesses (e.g., MuzeroUe et al. 2006, Sicilia-Aguilar et al. 2006a) This definition 
has been relaxed by some authors (e.g.. Brown et al. 2007, Cieza et al. 2008) to include 
objects with small, but still detectable, near-IR excesses. Transition disks have also been 
more broadly defined in terms of a significant decrement relative to the Taurus median 
Spectral Energy Distribution (SED) at any or all wavelengths (e.g., Najita et al. 2007). 
This broad definition is the closest one to the criteria we adopt to select our sample (see 
§ [2]). Even though, according to our definition, many transition disks have inner disks 
with non-zero opacity, we refer to their regions of low opacity as the "inner opacity hole." 
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Recent submillimeter studies provide dramatic support for the presence of inner opacity 
holes inferred from the SED modehng of transition disks. High resolution submillimeter 
continuum images of objects such as LkHa 330 (Brown et al. 2008) and GM Tau (Hughes 
et al. 2009) show sharp inner holes tens of AU in radius and lend confidence to the standard 
interpretation of transition disk SEDs. 

One of the most intriguing results from transition disk studies has been the great 
diversity of SED morphologies revealed by the Spitzer Space Telescope. In an attempt 
to describe distinct classes of transition disks, several names have recently emerged in 
the literature (Evans et al. 2009). These new names include: anemic disks, flat disks, 
or homologously depleted disks to describe objects whose observed SEDs are significantly 
below the median of the CTTS population at all IR wavelengths (Lada et al. 2006, Currie et 
al. 2009). These new names also include cold disks, objects with little or no near-IR excesses 
whose SEDs raise very steeply in the mid-IR (Brown et al. 2007), and pre-transitional disks, 
disks with evidence for an optically-thin gap separating optically-thick inner and outer disk 
components (Espaillat et al. 2008). 

Studying the diverse population of transition disks is key for understanding circumstellar 
disk evolution as much of the diversity of their SED morphologies is likely to arise from 
different physical processes dominating the evolution of different disks. Disk evolution 
processes include: viscous accretion (Hartmann et al. 1998), photoevaporation (Alexander 
et al., 2006), the magnetorotational instability (Chiang & Murray-Clay, 2007), grain growth 
and dust settling (Dominik & DuUemond, 2008), planet formation (Lissauer, 2003; Boss et 
al. 2000), and dynamical interactions between the disk and stellar or substellar companions 
(Artymowicz & Lubow, 1994). Understanding the relative importance of these physical 
processes in disk evolution and their connection to the different classes of transition disks is 
currently one of the main challenges of the disk evolution field. Also, even though transition 
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disks are rare, it is likely that all circumstellar disks go through a short transition disk stage 
(as defined by their SEDs) at some point of their evolution. This is so because observations 
show that an IR excess at a given wavelength is always accompanied by an excess at longer 
wavelengths (the converse is not true as attested by the very existence of transition disks) . 
This implies that, unless some disks manage to lose the near-, mid- and far-IR excess at 
exactly the same time, the near-IR excess always dissipates before the mid-IR and far-IR 
excess do. Since no known process is expected to remove the circumtellar dust at all radii 
simultaneously, it is reasonable to conclude that transition disks represent a common (if not 
unavoidable) phase in the evolution of a circumstellar disk. 

At least four different mechanisms have been proposed to explain the opacity holes 
of transition disks: giant planet formation, grain growth, photoevaporation, and tidal 
truncation in close binaries. However, since all these processes can in principle result in 
similar IR SEDs, additional observational constraints are necessary to distinguish among 
them. As discussed by Najita et al. (2007), Cieza (2008), and Alexander (2008), disk 
masses, accretion rates, and multiplicity information are particularly useful to distinguish 
between the different mechanisms that are likely to produce the inner holes in transition 
disks. A vivid example of the need for these kind of data is the famous disk around CoKu 
Tau/4. Since its sharp inner hole was discovered by Spitzer (Forrest et al. 2004), its origin 
has been a matter of great debate. The hole was initially modeled to be carved by a giant 
planet (Quillen et al. 2004), but Najita et al. (2007) argued that the low mass and low 
accretion rate of the CoKu Tau/4 disk are more consistent with photoevaporation than with 
a planet formation scenario. More recently, Coku Tau/4 has has been shown to be a close 
binary star system with a 8 AU projected separation, rendering its disk a circumbinary one 
(Ireland & Kraus 2008). 

Since the number of well characterized transition disks is still in the tens, most studies 
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so far have focused on the modehng of individual objects such as TW Hydra, (Calvet et 
al. 2004), GM Aur and DM Tau (Calvet et al. 2005), LkHa 330 (Brown et al. 2008), 
and LkCa 15 (Espaillat et al. 2008). To date, few studies have discussed the properties 
of transition disks as a group (e.g., Najita et al. 2007; Cieza et al. 2008). These papers 
studied relatively small samples, suffer from different selection biases, and arrived, not 
too surprisingly, at very different conclusions. Najita et al. (2007) studied a sample of 12 
transition objects in Taurus and found that they have stellar accretion rates ~10 times 
lower and a median disk mass (~25 Mjup) that is ~4 times larger than the rest of the 
disks in Taurus. They argue that most of the transition disks in their sample are consistent 
with the planet formation scenario. The disk masses found by Najita et al. are in stark 
contrast to the results from the SMA study of 26 transition disks from Cieza et al. (2008). 
They observed mostly WTTSs disks and found that all of them have very low masses <l-3 
Mjup suggesting that their inner holes were more likely due to photoevaporation, instead 
of the formation of jovian planets. This discrepancy can probably be traced back to the 
different sample selection criteria as Najita et al. studied mostly CTTSs, while Cieza et al. 



studied mostly WTTSs (see § 5.1.1). A much larger and unbiased sample of transition disks 
is needed to quantify the importance of multiplicity, photoevaporation, grain growth, and 
planet formation on the evolution of circumstellar disks. 

This paper is the first part of a series from an ongoing project aiming to characterize 
over 100 Spitzer-selected transition disks located in nearby star-forming regions. Here we 
present millimeter wavelength photometry (from the SMA and the CSO), high-resolution 
optical spectroscopy (from the Clay, CFHT, and Du Pont telescopes), and Adaptive 
Optics near-infrared imaging (from the VLT) for a sample of 26 Spitzer-selected transition 
circumstellar disks located in the Ophiuchus molecular cloud. We use these new ground- 
based data to estimate the disk mass, accretion rate, and multiplicity for each object in 
our sample in order to investigate the mechanisms potentially responsible for their inner 
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opacity holes. The structure of this paper is as follows. Our sample selection criteria are 
presented in § 2, while our observations and data reduction techniques are described in 
§ 3. We present our results on disk masses, accretion rates, and multiplicity in § 4. In 
§ 5, we discuss the properties of our transition disk sample and compare them to those of 
non-transition objects. We also discuss the likely origins of the inner holes of individual 
targets and the implications of our results for disk evolution. Finally, a summary of our 
results and conclusions is presented in § 6. 



2. Sample Selection 

We drew our sample from the 297 Young Stellar Object Candidates (YSOc) in the 
Ophiuchus catalog of the Cores to Disks (Evans et al. 2003) Spitzer Legacy Project. For a 
description of the Cores to Disks data products, see Evans et al. (2007) In particular, we 
selected all the targets meeting the following criteria: 

a. Have Spitzer colors [3.6]-[4.5] < 0.25. These YSOc are objects with small or no near-IR 
excess (see Figure Q. The lack of a [3.6] - [4.5] color excess in our targets is inconsistent 
with an optically thick disk extending inward to the dust sublimation radius, and therefore 
indicates the presence of an inner opacity hole. The presence of this inner opacity hole is 
the defining feature we intend to capture in our sample. This feature is present in ~21% of 
the YSOc in Ophiuchus as our first criterion selects 63 of them (out of 297). 

b. Have Spitzer colors [3.6]-[24] > 1.5. We apply this criterion to ensure all our targets have 
very significant excesses. It removes the 10 YSOc with smallest 24 yum excess (i.e., 1.5 > 
[3.6]-[24] > 1.0) and leaves 53 targets. 



^available at the Infra-Red Science Archive 


http : //irsa . ipac . caltech . edu/data/SPITZER/C2D/ 


^available at 


http: //irsa. ipac . caltech 


. edu 


/data /SPITZER/C2D/doc 
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c. Are detected with signal to noise ratio > 7 in all 2MASS and IRAC wavelengths as well 
as at 24 fim, to ensure we only include targets with very reliable photometry. This criterion 
removes 9 objects and leaves 44 targets. 

d. Have < 11 mag , driven by the sensitivity of our near-IR adaptive optics observations 
and to ensure a negligible extragalactic contamination (Padgett et al. 2008). This criterion 
removes 6 objects and leaves 38 targets. 

e. Are brighter than R = 18 mag according to the USNO-Bl catalog (Monet et al. 2003), 
driven by the sensitivity of our optical spectroscopy observations. This criterion removes 4 
objects and leaves a final target list of 34 YSOc. 

The first two selection criteria ( [3.6]-[4.5] < 0.25 and [3.6]-[24] > 1.5 ) effectively 
become our working definition for a transition disk. These criteria are fairly inclusive and 
encompass most of the transition disk definitions discussed in § [T] as they select disks with a 
significant flux decrement relative to "full disks" in the near-IR or at all wavelengths. Many 



of our targets have IR excesses that only become evident at 24 /xm (see § 5.1.3). To check 
the reality of these excesses, we have visually inspected the cutouts of the 24 /im Ophiuchus 
mosaic created by the Cores to Disks project]^ to confirm they show bona fide detections 
of our targets. We have also verified that all the 24 /xm images of our targets have been 
assigned an "Image Type" = 1 in the Cores to Disks catalogs, corresponding to objects 
that are well fitted by a point source profile. As a final check, we have verified that all our 
targets have "well behaved" SEDs that are consistent with reddened stellar photospheres 
shortward of 4.5 /im and IR excesses from a disk at longer wavelengths. We have observed 



all 34 YSOc in our target list. However, as discussed in § 4.1.2 this list includes one likely 



classical Be star and 7 likely Asymptotic Giant Branch stars. The remaining 26 targets 



available at 



http : //irsa. ipac . caltech.edu/data/SPITZER/C2D/ index_cut outs .html 
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are bona fide PMS stars with circumstellar disks and constitute our science sample. The 
Spitzer and ahernative names names, 2MASS and Spitzer fluxes, and the USNO-Bl R-band 
magnitudes for all our 34 targets are listed in Table 1. 



3. Observations 

3.1. Millimeter Wavelength Photometry 

Two of our 34 targets, #14 and 17, have already been detected at millimeter wavelength 
(Andrews & Williams, 2007 ), while stringent upper limits exist for 3 others, #12, 13, 
and 27 (Cieza et al. 2008). We have observed 24 of the 29 remaining objects with the 
Submillimeter Array (SMA; Ho et al. 2004), and 5 of them with Bolocam at the Caltech 



Submillimeter Observatory (CSO). In § 4.2, we use the millimeter wavelength photometry 



to constrain the masses of our transition disks. 



3.1.1. Submillimeter Array Observations 

Millimeter interferometric observations of 24 of our targets were conducted in service 
mode with the SMA, on Mauna Kea, Hawaii, during the Spring and Summer of 2009 (April 
Qth tj-^rough July IG*'') in the compact-north conflguration and with the 230 GHz/1300 fim 
receiver. Both the upper and lower sideband data were used, resulting in a total bandwidth 
of 4 GHz. 

Typical zenith opacities during our observations were T225ghz ~0. 15-0.25. For each 
target, the observations cycled between the target and two gain calibrators, 1625-254 
and 1626-298, with 20-30 minutes on target and 7.5 minutes on each calibrator. The 
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raw visibility data were calibrated with the MIR reduction packagqj The passband was 
flattened using ~1 hour scans of 3c454.3 and the solutions for the antenna-based complex 
gains were obtained using the primary calibrator 1625-254. These gains, applied to our 
secondary calibrator, 1626-298, served as a consistency check for the solutions. 

The absolute flux scale was determined through observations of either Callisto or 
Ceres and is estimated to be accurate to 15%. The flux densities of detected sources were 
measured by fltting a point source model to the visibility data, while upper limits were 
derived from the rms of the visibility amplitudes. The rms noise of our SMA observations 
range from 1.5 to 5 mJy per beam. We detected, at the 3-ct level or better, 5 of our 24 
SMA targets: #3, 15, 18, 21, and 32. The 1.3 mm fluxes (and 3-sigma upper limits) for our 
entire SMA sample are listed in Table 2. 



3.1.2. CSO-Bolocam Observations 

Millimeter wavelength observations of 5 of our targets were made with Bolocam[^ at the 
CSO on Mauna Kea, Hawaii, during June 25*'^-30*'^, 2009. The observations were performed 
in the 1.1 mm mode, which has a bandwidth of 45 GHz centered at 268 GHz. Our sources 
were scanned using a Lissajous pattern providing 10 min of integration time per scan. 
Between 16 and 40 scans per source were obtained. The weather was clear for the run, 
with r225GHz ranging from 0.05 to 0.1. Several quasars close to the science flelds were used 
as poniting calibrators, while Uranus and Neptune were used as flux calibrators. The data 
were reduced and the maps from individual scans were coadded using the Bolocam analysis 



"^available at 


http : //cf a-www . harvard . edu/~cqi/mircook . html 
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http: //www. CSO . caltech. edu/bolocam/ 
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pipeline |^ which consists of a series of modular IDL routines. None of the 5 targets were 
detected by Bolocam in the coadded maps, which have rms noises ranging from 3 to 6 mJy 
per beam. The 3-0" upper limits for the 1.1 mm fluxes of our Bolocam targets are listed in 
Table 2. 



3.2. Optical Spectroscopy 

We obtained Echelle spectroscopy (resolution > 20,000) for our entire sample using 3 
different telescopes. Clay, CFHT, and Du Pont. All the spectra include the Ha line, which 



we use to derive accretion rates (see § 4.3). 



3.2.1. Clay-Mike Observations 

We observed 14 of our targets with the Magellan Inamori Kyocera Echelle (MIKE) 
spectrograph on the 6.5-m Clay telescope at Las Campanas Observatory, Chile. The 
observations were performed in visitor mode on April 27*^^ and 28*^^, 2009. We used the red 
arm of the spectrograph and a 1" slit to obtain complete optical spectra between 4900 and 
9500 A at a resolution of 22,000. This resolution corresponds to ~0.3 A at the location of 
the Ha line and to a velocity dispersion of ~14 km/s. Since the CCD of MlKE's red arm 
has a pixel scale of 0.05 A/pixel, we binned the detector by a factor of 3 in the dispersion 
direction and a factor of 2 in the spatial direction in order to reduce the readout time and 
noise. The R-band magnitudes of our MIKE targets range from 15.5 to 18. For each object, 
we obtained a set of 3 or 4 spectra, with exposure times ranging from 3 to 10 minutes each, 
depending of the brightness of the targets. The data were reduced using the standard IRAF 

^ http: //www. cso . caltech. edu/bolocam/AnalysisSof tware .html 
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packages IMRED:CDDRED and ECHELLE:DOECSLIT. 

3.2.2. CFHT-Espadons Observations 

Twelve of our targets were observed with the ESPaDonS Echelle spectrograph on 
the 3.5-m Canada-France-Hawaii Telescope (CFHT) at Mauna Kea Observatory. The 
observations were performed in service mode during the last ESPaDonS observing run 
of the 2009A semester (June 30*'^-July IS*'^). The spectra were obtained in the standard 
"star+sky" mode, which delivers the complete optical spectra between 3500 A and 10500 A 
at a resolution of 68,000, or 4.4 km/s. The R-band magnitudes of our ESPaDonS targets 
range between 7 and 15. For each object, we obtained a set of 3 spectra with exposures 
times ranging from 2.5 to 10 minutes each, depending on the brightness of the target. The 
data were reduced through the standard CFHT pipeline Upena, which is based on the 
reduction package Libre-ESpRIlQ 

3.2.3. Du Pont-Echelle Observations 

We observed 8 of our targets with the Echelle Spectrograph on the 2.5-m Irenee 
du Pont telescope at las Campans Observatory. The observations were performed in 
visitor mode between May 14*^ and May 16*^, 2009. We used a 1" slit to obtain spectra 
between 4000 and 9000 A with a resolution of 45,000, corresponding to 0.14 A in the red. 
However, since the Spectrograph's CCD has a pixel scale of ~0.1 A/pixel, the true two-pixel 
resolution corresponds to ~32,000, or ~9.4 km/s in the red. The R-band magnitudes of 
our Du Pont targets range between 15.0 and 16.5. For each object, we obtained a set of 

^ http : //www. cfht .hawaii . edu/Instruments/Spectroscopy/Espadons/Espadons_esprit .html 



3 to 4 spectra with exposures times ranging from 10 to 15 minutes each, depending on 
the brightness of the target. The data reduction was performed using the standard IRAF 
packages IMRED:CDDRED and TWODSPEC:APEXTRACT. 



3.3. Adaptive Optics Imaging 



High spatial resolution near-IR observations of our entire sample were obtained with 
the Nasmyth Adaptive Optics Systems (NAOS) and the CONICA camera at the 8.2 m 
telescope Yepun, which is part the European Southern Observatory's (ESO) Very Large 
Telescope (VLT) in Cerro Paranal, Chile. The data were acquired in service mode during 
the ESO's observing period 84 (April I''* through September 30*'^, 2009). 

To take advantage of the near-IR brightness of our targets, we used the infrared 
wavefront sensor and the N90C10 dichroic to direct 90% of the near-IR light to adaptive 
optics systems and 10% of the light to the science camera. We used the S13 camera 
(13.3 mas/pixel and 14x14" field of view) and the Double RdRstRd readout mode. The 
observations were performed through the and J-band filters at 5 dithered positions per 
filter. The total exposures times ranged from 1 to 50 s for the K^-band observations and 
from 2 to 200 s for the J-band observations, depending on the brightness of the target. The 
data were reduced using the Jitter software, which is part of ESO's data reduction package 




In § 4.4, we use these Adaptive Optics (AO) data to constraint the multiplicity of 



our targets. 



http : //www . eso . org/pro j ects/ aot /eclipse/ 
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4. Results 

4.1. Stellair properties 

Before discussing the circumstellar properties of our targets, which is the main focus 
of our paper, we investigate their stellar properties. In this section, we derive their spectral 
types and identify any background objects that might be contaminating our sample of PMS 
stars with IR excesses. 

4.1.1. Spectral Types 

We estimate the spectral types of our targets by comparing temperature sensitive 
features in our echelle spectra to those in templates from stellar hbraries. We use the 
libraries presented by Soubiran et al. (1998) and Montes (1998). The former has a spectral 
resolution of 42,000 and covers the entire 4500 - 6800 A spectral range. The Latter has 
a resolution of 12,000 and covers the 4000 - 9000 A spectral range with some gaps in the 
coverage. Before performing the comparison, we normalize all the spectra and take the 
template and target to a common resolution. 

Most of our sources are M-type stars, for which we assign spectral types based on the 
strength of the TiO bands centered around 6300, 6700, 7150, and 7800 A. We classify G-K 
stars based on the ratio of the V I (at 6199 A) to Fe I (6200 A) hne (Padgett, 1996) and/or 
on the strength of the Ca I ( 6112 A) and Na I (5890 and 5896 A) absorption hues (Montes 
et al. 1999; Wilking et al. 2005). There is also a B-type star in our sample, which we 
identify and type by the width of the underlying Ha absorption line (which is much wider 
than its emission hne) and by the strength of the Paschen 16, 15, 14, and 13 hues. The 
spectral types so derived are listed in Table 2. We estimate the typical uncertainties in our 
spectral classification to be 1 spectral subclass for M-type stars and 2 spectral subclasses 
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for K and earlier type stars. 

4- 1-2. Pre-main-sequence stars identification 

Background objects are known to contaminate samples of Spitzer-selected YSO 
candidates (Harvey et al. 2007, Oliveira et al. 2009). At low flux levels, background galaxies 
are the main source of contamination. However, the optical and near-IR flux cuts we have 
implemented as part of our sample selection criteria seem to have very efficiently removed 
any extragalactic sources that could remain in the Cores to Disks catalog of Ophiuchus 
YSO candidates. At the bright end of the flux distribution. Asymptotic Giant Branch 
(AGB) stars and classical Be stars are the main source of contamination. AGB stars are 
surrounded by shells of dust and thus have small, but detectable, IR excesses. The extreme 
luminosites (~10'^ Lq) of AGB stars imply that they can pass our optical and near-lR 
flux cuts even if they are located several kpc away. Classical Be stars are surrounded by a 
gaseous circumstellar disk that is not related to the primordial accretion stage but to the 
rapid rotation of the object (Porter & Rivinius, 2003). Classical Be stars exhibit both IR 
excess (from free-free emission) and Ho; emission (from the recombination of the ionized 
hydrogen in the disk) and thus can easily be confused with Herbig Ae/Be stars (early-type 
PMS stars). 

There is only one B-type star in our target list, object # 2. This object is located 
close to the eastern edge of the Cores to Disks maps of Ophiuchus, shows very little 24 //m 
excess, and is not detected at 70 fim. Its weak IR excess is more consistent with the free-free 
emission of a classical Be star than with the thermal IR excess produced by circumstellar 
dust around a Herbig Ae/Be star. Since there is no clear evidence that object # 2 is in fact 
a pre-main-sequence star, we do not include it in our sample of transition disks. 
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G and later type PMS stars can be distinguished from AGB stars by the presence of 
emission hnes associated with chromospheric activity and/or accretion. Ha (6562 A) and 
the Ca II infrared triplet (8498, 8542, 8662 A) are the most conspicuous of such lines. The 
Li I 6707 A absorption line is also a very good indicator of stellar youth in mid-K to M-type 
stars because Li is burned very efficiently in the convective interiors of low-mass stars and is 
depleted soon after these objects arrive on the main-sequence. In Table 2, we tabulate the 
velocity dispersion of the Ha emission lines, the equivalent widths of the detected Li I 6707 
A absorption lines, and the identification of the Ca II infrared triplets. After removing the 
Be star from consideration, we detect Ha emission in 25 of our targets. All of them exhibit 
Li I absorption and/or Ca II emission and can therefore be considered bona fide PMS stars. 
Six of the targets exhibit Ha absorption. All of them are M-type stars with small 24 fim 
excesses and no evidence for Li I absorption nor Ca II emission. Since these targets are 
most likely background AGB stars, we do not include them in our sample of transition 
disks. Two of our targets, objects # 12 and 22, show no evidence for Ha emission or 
absorption. Object #12 is the well studied K-type PMS star Do Ar 21, and we thus include 
it in our transition disks sample. Object 7^ 22 is a M6 star with small 24 fim excesses and 
no evidence for Li I absorption nor Ca II emission and is also a likely AGB star. 

To sum up, after removing one classical Be star (object # 2) and seven likely AGB stars 
(objects # 4, 6, 8, 10, 22, 33, 34), we are left with 26 bona-fide PMS stars. These 26 objects 
constitute our sample of transition disks. As shown in Figure [T| all the contaminating 
objects have [3.6]-[24] < 2.5, while 24 of the 26 PMS stars have [3.6]-[24] > 2.5. This result 
underscores the need for spectroscopic diagnostics to firmly establish the PMS sequence 
status of Spitzer-selected YSO candidates, especially those with small IR excesses. 



4.2. Disk Masses 



As shown by Andrews & Williams (2005,2007), disk masses obtained from modeling the 
IR and (sub)-mm SEDs of circumstellar disks are well described by a simple linear relation 
between (sub)-mm flux and disk mass. From the ratios of model-derived disk masses to 
observed mm fluxes presented by Andrews & Williams (2005) for 33 Taurus stars, Cieza et 
al. (2008) obtained the following relation, which we adopt to estimate the disk masses of 
our transition disks: 



MnisK = 1.7 X io-i[( -^^(l-3^H ) ^ (^)2]M,^p (1) 

mJy 140pc 

Based on the standard deviation in the ratios of the model-derived masses to 
observed mm fluxes, the above relation gives disk masses that are within a factor of ~2 
of model-derived values; nevertheless, larger systematic errors can not be ruled out (see 
Andrews & Williams, 2007). In particular, the models from Andrews & Williams (2005, 
2007) assume an opacity as a function of frequency of the form oc u and a normalization 
of Kq = 0.1 gr/cm^ at 1000 GHz. This opacity implicitly assumes a gas to dust mass ratio 
of 100. Both the opacity function and the gas to dust mass ratio are uncertain and expected 
to change due to disk evolution processes such as grain growth and photoevaporation. 
Detailed modeling and a additional observational constraints on the grain size distributions 
and the gas to dust ratios will be needed to derive more accurate disk masses for each 
individual transition disk. 

Since in the (sub)mm regime disk fluxes behave as F^, oc i/2±o.5 (^Andrews & Williams, 
2005), our targets are expected to be brighter (by a factor of ~1.4) at 1.1 mm than they 
are at 1.3 mm. We thus modify Equation [T] accordingly to derive upper limits for the disk 
masses of the objects observed at 1.1 mm with the CSO. The disk masses (and 3-a upper 
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limits) for our sample, derived adopting a distance to Opliiuclius of 125 pc (Loinard et 
al. 2008), are listed in Table 3. For two objects, #12 and 27, we adopt the disk mass 
upper limits from Cieza et al. (2008), which were derived from 850 fim observations. The 
vast majority of our transition disks have estimated disk masses lower than ~2.0 Mjup- 
However, two of them have disk masses typical of CTTSs (~3-5 Mjup) and two others have 
significantly more massive disks (~10-40 Mjup). 



4.3. Accretion Rates 

Most PMS stars show Ha emission, either from chromospheric activity or accretion. 
Non-accreting objects show narrow (< 200 km/s) and symmetric line profiles of chromspheric 
origin, while accreting objects present broad (> 300 km/s) and asymmetric profiles produced 
by large- velocity magnetospheric accretion columns. The velocity dispersions (at 10 % 
of the peak intensity), AV, of the Ha emission lines of our transition disks are listed in 
Table 2. The boundary between accreting and non-accreting objects has been empirically 
placed by different studies at AV between 200 km/s (Jayawardana et al. 2003) and 270 
km/s (White & Basri, 2003). Since only one of our objects, source # 25, has AV^ in the 
200-300 km/s range, most accreting and non-accreting objects are clearly separated in our 
sample. Source # 25 has AV^ ~230 km/s and a very noisy spectrum. We consider it to 
be non-accreting because it has a very low fractional disk luminosity (L/j/L^, < 10^'^, see 



5.1.3). 



The continuum-subtracted Ha profiles for all the 17 accreting transition disks are 
shown in Figure [2| while those for the 8 non-accreting disks where the Ha line was detected 
are shown in Figure |3| For accreting objects, the velocity dispersion of the Ha line correlates 
well with accretion rates derived from detailed models of the magnetospheric accretion 
process. We therefore estimate the accretion rates of our targets from the width of the 
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Ha line measured at 10% of its peak intensity, adopting the relation given by Natta et al. 
(2004): 

Log{Macc{MQ/yr)) = -12.89(±0.3) + 9.7(±0.7) x 10-^AV{km/s) (2) 

This relation is valid for 600 km/s > AV > 200 km/s (corresponding to 10^'^ Mq/jt 
> Mace 10"^^ Mq/jt) and can be applied to objects with a range of stellar (and sub-stellar) 
masses. The broadening of the the Ha line is 1 to 2 orders of magnitude more sensitive 
to the presence of low accretion rates than other accretion indicators such as U-band 
excess and continuum veiling measurements (Muzerolle et al. 2003, Sicilia-Aguilar et al. 
(2006b), and is thus particularly useful to distinguish weakly accreting from non-accreting 
objects. However, as discussed by Nguyen et al. (2009), the 10% width measurements are 
also dependent on the line profile, rendering the 10% Ha velocity width a relatively poor 
quantitive accretion indicator, specially at high accretion rates (Muzerolle et al. 1998). 

For the 9 objects we consider to be non-accreting, we adopt a mass accretion upper 
limit of lO^^^M0/yr, corresponding to AV = 200 km/s. The so derived accretion rates 
(and upper limits) for our sample of transition disks are listed in Table 3. Given the large 
uncertainties associated with Equation 2 and the intrinsic variability of accretion in PMS 
stars, these accretion rates should be considered order-of-magnitude estimates. 

4.4. Stellar companions 

Seven binary systems were identified in our sample by visual inspection of the VLT-AO 
images: targets #5, 9, 15, 18, 23, 24, and 27 (Figure |4]). The separation and flux ratios 
of these systems range from 0.19" to 1.68" and 1.0 to 12, respectively. All these systems 
were well resolved in both our J- and K5-band images, which have typical FWHM values 
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of 0.06"-0.07". Targets #9, 15, 18, and 27 are previously known binaries (Ratzka et al. 
2005; Close et al. 2007), while targets #5, 23, and 24 are newly identified multiple systems. 
For each one of these binary systems, we searched for additional tight companions by 
comparing each other's PSFs. The PSF pairs were virtually identical in all cases, except for 
target #24. The north-west component of this target has a perfectly round PSF, while the 
south-east component, 0.84" away, is clearly elongated. Since variations in the PSF shape 
are not expected within such small angular distances and this behavior is seen in both the 
J- and Kg-band images, we conclude that target #24 is a triple system. We use the PSF 
of the wide component (object #24- A) to model the image of the tight components (object 
#24-B/C). We find that the image of the tight components is well reproduced by two 
equal-brightness objects 4 pixels apart (corresponding to ~0.05" or 6.6 AU) at a position 
angle of ~30 deg. 

We have also searched in the literature for additional companions in our sample 
that our VLT observations could have missed. In addition to the seven multiple systems 
discussed above, we find that source # 12 (Do Ar 21) is a binary system with a projected 
separation of just 0.005" (corresponding to 0.62 AU), recently identified by the Very Long 
Baseline Array (Loinard et al. 2008). We also find that source # 27 is in fact a triple 
system. The "primary" star in the VLT observations is itself a spectroscopic binary with 
a 35.9 d period and an estimated separation of 0.27 AU (Mathieu et al. 1994). Additional 
multiphcity constraints exists for 2 other of our targets. Target #21 has been observed 
with the lunar occultation technique (Simon et al. 1995) and a similar-brightess companion 
can be ruled out down to ~1 AU. Similarly, target #13 has been found to have a constant 
radial velocity ( cr^ei < 0.25 km/s) during a 3 yr observing campaign (Prato, 2007). 

For the apparently single stars, we estimated the detection limits at 0.1" separation 
from the 5-a noise of PSF-subtracted images. For these objects, however, there are no 
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comparison PSF's available to perform the subtraction. Instead, we subtract a PSF 
constructed by azimuthally smoothing the image of the target itself, as follows. For each 
pixel in the image, the separation from the target's centroid is calculated, with sub-pixel 
accuracy. The median intensity at that separation, but within an arc of 30 pixels in 
length, is then subtracted from the target pixel. Thus, any large scale, radially symmetric 
structures arc removed. The detection limits at 0.1" separation, obtained as described 
above, range from 2.4 to 3.5 mag, with a median value of 3.1 mag (corresponding to a flux 
ratio of 17). The separations, positions angles, and flux ratios of the multiple systems in 
our sample are shown in Table 2, together with the flux ratio hmits of companions for the 
targets that appear to be single. 

5. Discussion 

5.1. Sample properties 

In this section, we investigate the properties of transition disks as a group of objects. 
We discuss the properties of our targets (disk masses, accretion rates, multiplicity, SED 
morphologies, and fractional disk luminosities) and compare them to those of non-transition 
disks in order to place our sample into a broader context of disk evolution. 

5.1.1. Disk Masses and Accretion Rates 

To compare the disk masses and accretion rates of our sample of transition disks to 
those of non-transition disks, we collected such data from the literature for Ophichus objects 
with Spitzer colors [3.6]-[4.5] > 0.25 and [3.6]-[24] > 1.5. These are objects that have both 
near-IR and mid-IR excesses and therefore are likely to have "full disks" , extending inward 
to the dust sublimation radius. We collected disk masses for non-transition disks from 
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Andrews & Williams (2007). This comparison sample is appropriate because our mm flux 
to disk mass relation (i.e., Equation [T]) has been calibrated using their disk models. The 
accretion rates for our comparison sample come from Natta et. (2006), and have been 
calculated from the luminosity of the Pa/3 line. This comparison sample is also appropriate 
because the relations between the luminosity of the Pa/3 line and accretion rate and between 
the Ha full width at 10% intensity and accretion rate (which we use for our transition disk 
sample) have both been calibrated by Natta et al. (2004) using the same detailed models of 
magnethospheric accretion. 

Figure |5] shows that transition disks tend to have much lower disk masses and accretion 
rates than non-transition objects. We also flnd a strong connection between the magnitude 
of the mid-IR excess and both disk mass and accretion rate: objects with little (< 4 mag) 
24 /im excess tend to have small disk masses and low accretion rates. Among our transition 
disk sample, we also note that all of the mm detections correspond to accreting objects and 
that even some of the strongest accretors have very low disk masses (see Figure |6]). In other 
words, massive disks are the most likely to accrete, but some low-mass disks can also be 
strong accretors. 

These results can help reconcile the apparently contradictory flndings of previous 
studies of transition disks. As discussed in § [T| Najita et al. (2007) studied a sample 
of 12 transition objects in Taurus and found that they have a median disk mass of 25 
Mj[/p, which is ~4 times larger than the rest of the disks in Taurus, while Cieza et al. 
(2008) found that the vast majority of their 26 transition disk targets had very small (< 2 
Mj(/p) disk masses. Our results now show that transition disks are a highly heterogeneous 
group of objects, whose "mean properties" are highly dependent on the details of the 
sample selection criteria. On the one hand, the sample studied by Cieza et al. (2008) was 
dominated by weak-line T Tauri stars. This explains the low disk masses they found, as 
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lack of accretion is a particularly good indicator of a low disk mass (See Figure [6]). One the 
other hand, Najita et al. (2006) drew their sample from the Spitzer spectroscopic survey 
of Taurus presented by Furlan et al. (2006), which in turn selected their sample based on 
mid-IR colors and fluxes from the Infrared Astronomical Satellite (IRAS). As a result, the 
Spitzer spectroscopic survey of Taurus was clearly biased towards the brightest objects in 
the mid-IR. This helps to explain the high disk masses obtained for the transition disks 
studied by Najita et al. (2006), as objects with strong mid-IR excesses tend to have higher 
disk masses (Figure [s]). 



5.1.2. Multiplicity 

As discussed in § [T| dynamical interactions between the disk and stellar companion is 
one of the main mechanisms proposed to account for the opacity holes of transition disks. 
However, the fraction of transition disks that are in fact close binaries still remains to 
be established. If binaries play a dominant role in transition disks, one would expect a 
higher incidence of close binaries in transition disks that in non-transition disks. In order 
to test this hypothesis, we compare the distribution of binary separations of our sample 
of transition disk to that of non-transition disk (and disk- less PMS stars). We drew our 
sample of non-transition disks and disk-less PMS stars from the compilation of Ophiuchus 
binaries presented by Cieza et al. (2009). 

As shown by Figure [7], disk-less PMS stars tend to have companions at smaller 
separations than stars with regular, non-transition disks. According to a two-sided 
Kolmogoro-Smirnov (KS) test, there is only a 0.5% probability that the distributions of 
binary separations of non-transition disks and disk-less stars have been drawn from the 
same parent population. As discussed by Cieza et al. (2009), this result can be understood 
in terms of the effect binaries have on circumstellar disks lifetimes via the tidal truncation 
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of the outer disk. For instance, a binary system with a 30 AU separation is expected 
to initially have individual disks that are ~10-15 AU in radius. Given that the viscous 
timescale is roughly proportional to the size of the disk, such small disks are likely to have 
very small accretion lifetimes, smaller than the age of the sample. 

For systems with binary separations much smaller than the size of a typical disk, 
the outer disk can survive in the form of a circumbinary disk, with a tidally truncated 
inner hole with a radius ~2x the orbital separation (Artymowicz & Lubow, 1994). Such 
systems could be classified as transition disks based on their SEDs. We do not see an 
increased incidence of binaries with separations in the ~8-20 AU range, the range where 
we are sensitive to companions that could in principle carve the inferred inner holes. In 
fact, the incidence of binaries in this separation range is lower for our transition disks 
than it is for the samples of both non-transition disks and disk-less stars. This suggests 
that stellar companions at 8-20 AU separations are not responsible for a large fraction of 
the transition disk population. Given the size of our transition disk sample, however, this 
result carries little statistical weight, and should be considered only a hint. According to 
a KS test, the distribution of binaries separations of transition objects is indistinguishable 
from that of both non-transition disks (p = 28%) and disk-less stars (p = 12%). On the 
other hand, the fact that 2 of our transition disk targets have previously known companions 
at sub-AU separations suggests that circumbinary disks around very tight systems (e.g., 
r < 1-5 AU) could indeed represent an important component of the transition disk 
population. Unfortunately, the presence of tight companions (r < 8 AU) remains completely 
unconstrained for most of our sample. A complementary radial velocity survey to find the 
tightest companions is highly desirable to firmly establish the fraction of transition disks 
that could be accounted for by close binaries. 
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5.1.3. SED morphologies and fractional disk luminosities 

In addition to the disk mass, accretion rate, and multiplicity, the SED morphology 
and fractional disk luminosity of a transition disk can provide important clues on the 
nature of the object. In this section, we quantify the SED morphologies and fractional disk 
luminosities of our transition objects in order to compare them to the values of objects with 
"full disks". 

The diversity of SED morphologies seen in transition disks can not be properly 
captured by the traditional classification scheme of young stellar objects (i.e., the class I 
through III system), which is based on the slope of the SED between 2 and 25 /xm (Lada 
1987). Instead, we quantify the SED "shape" of our targets adopting the two-parameter 
scheme introduced by Cieza et al. (2007), which is based on the longest wavelength at 
which the observed flux is dominated by the stellar photosphere, Xtum-off, and the slope 
of the IR excess, aexcess, computed as dlog(AF)/dlog(A) between Xtum-off and 24 /^m. The 
Xturn-off and (Xexcess valucs for our entire sample are listed in Table 3. 

The Xturn-off and aexcess parameters provide useful information on the structure of 
the disk. In particular, Xturn-off clearly correlates with the size of the inner hole as it 
depends on the temperature of the dust closest to the star. Wien's displacement law implies 
that Spitzer's 3.6, 4.5, 5.8, 8.0, and 24 fim bands best probe dust temperatures of roughly 
780K, 620K, 480K, and 350K, and 120K, respectively. For stars of solar luminosity, these 
temperatures correspond to circumstellar distances of <1 AU for the IRAC bands and 
~10 AU for the 24 /im MIPS band. Similarly, aexcess correlates well with the sharpness of 
the hole: a sharp inner hole, such as that of CoKu Tau/4, is characterized by a large and 
positive aexcess value, while a continuous disk that has undergone signiflcant grain growth 
and dust settling is characterized by a large negative aexcess value (DuUemond & Dominik, 
2004). However, the Xtum-off and aexcess parameters should be interpreted with caution. 
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Deriving accurate disk properties from a SED requires detailed modeling because other 
properties, such as stellar luminosity and disk inclination, can also affect the Xtum-off and 

dexcess ValuCS. 

The fractional disk luminosity, the ratio of the disk luminosity to the stellar luminosity 
(L/j/L*), is another important quantity that relates to the evolutionary status of a disk. 
On the one hand, typical primordial disks around CTTSs have L/)/L* ~ 0.1 as they have 
optically thick disks that intercept (and reemit in the IR) ~10% of the stellar radiation. 
On the other hand, debris disks show L/j/L^, values < 10~^ because they have optically thin 
disks that intercept and reprocess ~10~^-10~'^ of the star's light (e.g., Bryden et al. 2006). 

We estimate L^j/L* for our sample of transition disk^ using the following procedure. 
First, we estimate L^, as in Cieza et al (2007), by integrating over the broad-band colors 
given by Kenyon & Hartmann (1995) for stars of the same spectral type as each one of 
the stars in our sample. The integrated fluxes were normalized to the J-band magnitudes 
corrected for extinction, adopting Aj =1.53 x E(J-K5 ), where E(J-K5) is the observed 
color excess with respect to the expected photospheric color. Then, we estimated Ld 
by integrated the estimated disk fluxes, the observed fluxes (or upper limits) minus the 
contribution of stellar photosphere at each wavelength between Xtum-off and 1.3 mm (or 
850 fim for targets #12 and 27). The disk was assumed to contribute 30% of the flux 
at Xtiirn— off (e.g., a small but non- negligible amount, consistent with our definition of 

Xturn—off)- 

The near and mid-IR luminosities of our transition disks are well constrained because 
their SEDs are relatively well sampled at these wavelengths. However, their far-IR 

^for completeness, we also estimate analogous L/j/L* values for the likely Be and AGB 
stars in our sample. 
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luminosities remain more poorly constrained. Only 9 of our 34 targets have ( 5-a or better) 
detections at 70 fim listed in the Cores to Disks catalog. For the rest of the objects, we 
have obtained 5-0" upper limits, from the noise of the 70 fim images at the location of 
the targets, in order to fill the gap in their SEDs between 24 fim and the millimeter. In 
particular, we used 120" x 120" cut outs^ of the Ophichus mosaic created by the Cores 
to Disks project from the "filtered" Basic Calibrated Data. We calculated the noise as 
1-cr = 2xl.7(rmSsA,.j/)v^, where imssky is the flux rms (in mJy) of an annulus centered on 
the source with an inner and outer radius of 40 and 64", respectively. N is the number 
of (4") pixels in a circular aperture of 16" in radius, and 1.7 is the aperture correction 
appropriate for such sky annulus and aperture size The factor of 2 accounts for the fact 
that the images we use have been resampled to pixels that are half the linear size of the 



original pixels, artificially reducing the noise of the images (by y where N^es and l^orig 
are, respectively, the number of resampled and original pixels in the sky annulus). 

The LOG(L£)/L*) values for our transition disk sample, ranging from -4.1 to -1.5 
are listed in Table 3. As most of the luminosity of a disk extending inward to the dust 
sublimation temperature is emitted in the near-IR, L^/L^, is a very strong function of 
Xturn-off- the shorter the Xtum-off wavelength, the higher the fractional disk luminosity. 
For objects with Xtum-off < 8-0 /im, the 70 fim flux represent only a minor contribution 
to the total disk luminosity. On the other hand, objects with with Xtum-off = 8.0 fim 
(i.e., objects with significant excess at 24 fim only) have much lower L^i/L* values, and the 
70 fim emission becomes a much larger fraction of the total disk luminosity. As a result, 
the Ld/^i* values of objets with Xtum-off = 8.0 /im and no 70 /im detections should be 
considered upper limits. 



^°available atlhttp: //irsa. ipac . caltech.edu/data/SPITZER/C2D/ index_cutouts .html 

^^http: / / ssc.spitzer.caltech.edu/mips/apercorr/ 
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As shown in Figure [8| our transition disks span the range Xtum-off = 2.2 to 8.0 yum 
and a^xcess ~ -2.5 to 1, while typical CTTSs occupy a much more restricted region of the 
parameter space. Based on the median and quartile SEDs presented by Furlan et al (2006), 
50% of late-type (K5-M2) CTTSs have Xtum-off = 1-25 /xm, -1.0 > aexcess > -0.64, and 
fractional disk luminosities in the 0.07 to 0.15 range. As discussed by Cieza et al. (2005), 
CTTS are likely to have J-band excesses. However, since this excess is at the < 30-40% 
level, the observed J-band fluxes are still dominated by the stellar photospheres, satisfying 
our Xtum-off definition. Figure [8] also shows that, as expected, transition disks with 
different L/j/L* values occupy different regions of the aexcess vs. Xtum-off space. Objects 
with L/j/L* > 10^^ lay close to the CTTSs loci, while objects with L/j/L* < 10^'^ all have 
Xtum-off = 8.0 /im and aexcess ^0.0. The objects in the latter group are all non- accreting. 
Their very low fractional disk luminosities and lack of accretion suggest they are already in 



a debris disk stage (see § 5.2.1). 



5.2. The Origin of the Opacity Holes 

The four different mechanisms that have been proposed to explain the inner holes of 
transition disks can be distinguished when disk masses, accretion rates and multiplicity 
information are available (Najita et al., 2007; Cieza, 2008; Alexander, 2008). In this section, 
we compare the properties of our transition disk sample to those expected for objects 
whose opacity holes have been produced by these four mechanisms: photoevaporation, tidal 
truncation in close binaries, grain-growth, and giant planet formation. 
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5.2.1. Photoevaporation 

According to photoevaporation models (e.g., Alexander et al. 2006), extreme UV 
photons, originating in the stellar chromosphere, ionize and heat the circumstellar hydrogen. 
Beyond some critical radius, the thermal velocity of the ionized hydrogen exceeds its 
escape velocity and the material is lost as a wind. At early stages in the evolution of the 
disk, the accretion rate across the disk dominates over the evaporation rate, and the disk 
undergoes standard viscous evolution. Later on, as the accretion rate drops to the level of 
the photoevaporation rate, the outer disk is no longer able to resupply the inner disk with 
material, the inner disk drains on a viscous timescale, and an inner hole is formed. Once 
an inner hole has formed, the entire disk dissipates very rapidly from the inside out as a 
consequence of the direct irradiation of the disk's inner edge. Thus, transition disks created 
by photoevaporation are expected to have relatively low masses {M^usk < ^^jup) and 
negligible accretion (<10^^° MQyr~^). 

We find that 9 of our 26 transition disks have low disk mass (< 0.6-2.5 Vijup) and 
negligible accretion (M^cc < 10^^^ MQ/year) and are hence consistent with photoevaporation. 
The SEDs of these objects are shown in Figure |9} These objects appear all in the lower left 
corner of Figure |6j which shows the disk mass as a function of accretion rate for our sample. 
However, they occupy different regions of the aexcess vs Xtum-off parameter space. The 
majority of these objects have aexcess ^ -l? the typical value of CTTS with "full disks" (see 



Figure 10). These objects are therefore consistent with relatively flat, radially continuous 
disks. Two objects, targets # 5 and 19, have aexcess ^ 0, and are thus suggestive of sharp 
inner holes. Target #12 also has an steep SED, but between 24 and 70 fim, suggesting a 
larger inner hole Sharp inner holes could also arise from dynamical interactions between 



^^The mid-IR excess of target #12 (DoAr 21) comes from heated material that is at least 
100 AU away from the star and it has been suggested it may originate from a small-scale 
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the disk and large bodies within it. Photoevaporation and dynamical interactions are 
not mutually exclusive processes, as photoevaporation can also operate on a dynamically 
truncated disk as long as the disk mass and accretion rates are low enough. Coku Tau/4 
is perhaps the most extreme example of this scenario. It is a 8 AU separation, equal-mass, 
binary system (Ireland & Kraus, 2008) resulting on a aexcess value of ~2. Just like targets 
#5, 12, and 19, Coku Tau/4 lacks accretion and has a very low disk mass (Najita et al. 
2007), and is consequently also susceptible to photoevaporation. Our VLT-AO observations 
reveal no stellar companion within ~8 AU of targets ^^5, 12, or 19. However, target #12 is 
already known to be very tight binary (r < 1 AU), and the presence of stellar (or planetary 
mass) companions remains completely unconstrained for targets ^^5 and 19. 

Photoevaporation provides a natural mechanism for objects to transition from the 
primordial to the debris disk stage. Once accretion stops, photoevaporation is capable of 
removing the remaining gas in a very short (<< 1 Myr) timescale (Alexander et al. 2006). 
As the gas in the disk photoevaporates, it is likely to carry with it the smallest grains 
present in the disk. What is left represents the initial conditions of a debris disk: a gas 
poor disk with large grains, planetesmals and/or planets. The very short photoevaporation 
timescale implies that some of the systems shown in Figure |9] could already be in a debris 
disk stage. In fact, the 4 non-accreting objects with L^i/L* < 10^^ have properties that 
are indistinguishable from those of young debris disks like AU Mic and GJ 182 (Liu et al. 
2004). These 4 objects are # 12, 20, 23, and 25. These objects might have already formed 
planets or might be in a multiple system (like target #12 is). Either way, most of the 
circumstellar gas and dust has already been depleted, and what we are currently observing 
could be the initial architectures of different debris systems, now subjected to processes 

photodissociation region powered by DoAr 21 itself and not necessarily from a circumstellar 
disk (Jensen et al. 2009). 
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such as dynamical interactions, the Poynting-Robertson effect and radiation pressure. The 
rest of the non accreting objects, #5, 9, 13, 19, and 24, have L^j/L^ > 10~^ and could be in 
the process of being photoe vapor ated. 

5.2.2. Close binaries 

Circumbinary disks are expected to be tidally truncated and have inner holes with a 
radius ~2x the orbital separation (Artymowicz & Lubow, 1994). It is believed that most 
PMS stars are in multiple systems (Ratzka et al. 2005) with the same semi-major axis 
distribution as MS solar-type stars (e.g., a lognormal distribution centered at ~28 AU). As 
a result, ~30% of all Ophiuchus binaries are expected to have separations in the ~l-20 AU 
range. If the primodial disk has survived in such systems, it is likely to be in the form of a 
circumbinary disk with tidally truncated inner holes up to 40 AU in radius. These close 
binary systems could hence be classified as transition disks based on their IR SEDs. 

Four of our targets have companions close enough to suspect they could be responsible 
for the observed hole, namely sources # 12, 27, 24, and 23, in increasing order of projected 
separation. Source # 12 is a binary with a projected separation of ~0.6 AU and Xturn-off— 
8.0 /im, which implies the disk must be a circumbinary one. However, given the long 
Xturn-off wavelength and tightness of the binary system, it seems unlikely that the present 
size of the inner hole is directly connected to the presence of the stellar companion. Source 
# 27 is triple system with projected separations of ~0.3 AU and 41 AU and Xtum-off — 
4.5 //.m. The IR excess could originate at a circumbinary disk around the tight components 
of the system or at a circumstellar disk around the wide component. Target #24 is also a 
triple system with projected separations of ~7 and 105 AU and Xturn-off — 5.8 /im. The IR 
excess could also in principle originate at a circumbinary disk around the tight components 
of the system. However, a circumprimary disk around source # 24-A is an equally likely 
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possibility. Target # 23 is a ~0.19" (or 24 AU) separation binary with Xturn-off — 8-0 
/xm. This object has a small 24 /xm excess, consistent with Wien side of the emission from 
a cold outer disk. Objects # 12 and 23 have a low disk masses (< 1.7 Mjup), negligible 
accretion, and L/j/L* < 10~^. As discussed in the previous section, they are likely to be in 
the debris disk stage, regardless of the origin of the inner hole. After target #23, the next 
tighter binary is source #5, with a projected separation of 0.54" (or 68 AU). Since this 
object has Xturn-off — 5.8 iim, suggesting the presence of dust at <1 AU distances, it is 
very unlikely that the presence of its inner hole is related to the observed companion. Since 
VLT-AO observations are clearly not sensitive to companions with separations < 8 AU, 
a complementary radial velocity survey to find the tightest companions would be highly 
desirable to identify additional circumbinary disks within our sample. 

5.2.3. Grain growth 

Once primordial sub-micron dust grains grow into somewhat larger bodies (r dust 
^ Xsteiiar-photons) , Hiost of the solid mass ceases to interact with the stellar radiation, 
and the opacity function decreases dramatically. Grain growth is a strong function of 
radius; it is more efficient in the inner regions where the surface density is higher and the 
dynamical timescales are shorter, and hence can also produce opacity holes. Idealized dust 
coagulation models (i.e., ignoring fragmentation and radial drift) predict extremely efficient 
grain growth (DuUemond & Dominik, 2005) resulting in the depletion of all small grains 
in timescales of the order of 10^ yrs, which is clearly inconsistent with the observational 
constraints on the lifetime of circusmtellar dust. In reality, however, the persistence of 
small opacity-bearing grains depends on a complex balance between dust coagulation and 
fragmentation (Dominik & DuUemond, 2008). 

Since grain growth affects only the dust and operates preferentially at smaller radii, a 
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disk evolution dominated by grain growth is expected to result in an actively accreting disk 
with reduced opacity in the inner regions. All the accreting disks in our sample, whose SEDs 



are shown in Figure 11, are thus in principle consistent with grain growth. Objects with 
very strong accretion (i.e., > 10~^ MqJt^'^) are especially good candidates for grain growth 
dominated disks as an opacity hole is likely to trigger the onset of the magneto-rotational 
instability and exacerbate accretion (Chiang & Murray-Clay, 2007). Target 7^ 14, which is 
the extensively studied system DoAr 25, is a prime example of this scenario. It has, by far, 
the highest disk mass in our sample (~38 Mj^p) and one of the highest accretion rates 
(~10^'''^ M0yr~^). This object has recently been imaged at high spatial resolution with the 
Submillimeter Array, and its SED and visibilities have been successfully reproduced with 
a simple model incorporating significant grain growth in its inner regions (Andrews et al. 
2008). 

Grain growth is considered the first step toward planet formation. Unfortunately, the 
observational effects of grains growing into terrestrial planets are no different from those 
of growing into meter size objects. As a result, the observations presented herein place no 
constraints on how far along the planet formation process currently is, unless the planet 
becomes a giant planet, massive enough to dynamically open a gap in the disk. 



5.2.4- Giant planet formation 

Since theoretical models of the dynamical interactions of forming giant planets with the 
disk (Lin & Papaloizous 1979, Artymowicz & Lubow 1994) predict the formation of inner 
holes and gaps, planet formation quickly became one of the most exciting explanations 
proposed for the inner holes of transition disks (Calvet et al. 2002, 2005; Quillen et al. 
2004; D'Alessio et al. 2005; Brown et. 2008). A planet massive enough to open a gap in 
the disk (M >0.1-0.5 Mjup), is expected to divert most of the material accreting from the 
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outer disk onto itself. As a result, in the presence of a Jupiter mass planet, the accretion 
onto the star is expected to be reduced by a factor of 4 to 10 with respect to the accretion 
across the outer disk (Lubow & D'Angelo, 2006). 

Four of the accreting objects in our sample, targets # 11, 21, 31, and 32 have aexcess 
> 0, suggesting the presence of a sharp inner hole, and are thus excellent candidates for 
ongoing giant planet formation. Source #32 is a particularly good candidate to be currently 
forming a massive giant planet, as it has a high disk mass (~11 Mjup) and small accretion 
rate (lO'^-^ MqJt-^). Sources #11, 21, and 31 (M^ < 1.5 Mjup; M„,e -lO'^-^-lO-^-^ 
M0yr~^), could have formed a massive giant planet in the recent past, as in their case 
most of the disk mass has already been depleted. Since target #21 has been observed with 
the lunar occultation technique, an equal-mass binary system can be ruled out for this 
system down to ~1 AU (Simon et al. 1995). Also, it has been proposed that companions 
with masses of the order of 10 Mjup completely isolate the inner disk from the outer disk 
and halt the accretion onto the star (Lubow et al. 1999). If that is the case, accretion 
itself could be evidence against the presence of a close (sub)stellar companion. That said, 
because a few accreting close binary systems are already known, such as DQ Tau (Carr 
et al. 2001) and CS Cha (Espaillat et al. 2007), it is clear that accretion does not rule 
out completely the presence of tight stellar companions. In addition to the mass of the 
companion, the eccentricity of the orbit, the viscosity and scale height of the disk are also 
important factors in determining whether or not gap-crossing streams can exist and allow 
the accretion onto the star to continue (Artymowicz &Lubow 1996). 



5.2.5. Disk Classification 



Based on the discussion above, we divide our transition disk sample into the following 
"disk types" : 
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a) 13 grain growth-dominated disks (accreting objects with a < 0). 

b) 4 giant planet forming disks (accreting objects with a > 0). 

c) 5 photoevaporating disks (non-accreting objects with disk mass < 2.5 Mjup, but L/^/L* 
> 10-=^) 

d) 4 debris disks (non-accreting objects with disk mass < 2.5 Mjup and L^/L* < 10~^) 

e) 4 circumbinary disks (a binary tight enough to accommodate both components within 
the inferred inner hole). 

The total number of objects listed add up to 30 instead of 26 because 4 objects fall 
into two categories. Sources # 12 and 23 are considered both debris disks and circumbinary 
disks. Object #24 has been classified as both a circumbinary disk and a photoevaporating 
disk, while object # 27 has been classified as both a circumbinary disk and a grain growth 
dominated disk. These last two objects are triple systems and their classification depends 
on whether the IR excess is associated with the tight components of the systems or with 
the wide components. 

The "disk types" for our targets are hsted in Table 3. All of our objects should be 
considered to be candidates for the categories listed above. The current classification 
represents our best guess given the available data. Only followup observations and detailed 
modeling will firmly establish the true nature of each object. An important caveat of this 
classification is the lack of constraints for most of the sample on stellar companions within 
< 8 AU, a range where ~30% of all stellar companions are expected to lay (Duquennoy 
& Mayor, 1991; Ratzka et al. 2005). Future radial velocity observations are very likely to 
increase the number of objects in the circumbinary disk category. 

While it is tempting to speculate that the "disk types" a through d represent an 
evolutionary sequence, starting with a "full disk" (i.e., a "full disk" that undergoes grain 
growth followed by giant planet formation, followed by photoevaporation leading to a debris 
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or disk-less stage) it is clear that not all disks follow this sequence. For instance, many of 
the grain growth dominated disks have disk masses that are 10 times smaller than some 
of the giant planet forming disks (e.g., targets #1, 15, 16, and 26 vs. #32). These former 
systems have very little mass left in the disk (< 1 Mjup) and might never form giant 
planets. It it also clear that not all disks become transition disks in the same way. Some 
objects become transition disks (as defined by their SEDs) while the total mass of the 
disk is very large (e.g., targets #32 and 14, ~ 10-40 Mjup), but other objects retain 
"full disks" even when the disk mass is relatively small (M^ ~2 Mjup), such as ROXRl 
29 (Cieza et al. 2008). Since very little dust is needed to keep a disk optically thick at 
near and mid-IR wavelengths, it is expected that some objects will only become transition 
disks when their disk masses and accretion rates are low enough to become susceptible to 
photoevaporation (i.e., they will evolve directly from a "full disk" to a photoevaporating 
disk followed by a debris or disk-less stage). 



6. Summary and Conclusions 

We have obtained millimeter wavelength photometry, high-resolution optical 
spectroscopy, and Adaptive Optics near-infrared imaging for a sample of 34 Spitzer-selected 
YSO candidates located in the Ophiuchus molecular cloud. All our targets have SEDs 
consistent with circumstellar disks with inner opacity holes (i.e., transition disks). After 
removing one likely classical Be star and seven likely AGB stars we were left with a sample 
of 26 transition disks. We have used these data to estimate the disk mass, accretion rate, 
and multiplicity of each transition disk in our sample in order to investigate the mechanisms 
potentially responsible for their inner opacity holes: dynamical interaction with a stellar 
companion, photoevaporation, grain growth, and giant planet formation. 

We find that transition disks exhibit a wide range of masses, accretion rates, and 
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SED morphologies. They clearly represent a heterogeneous group of objects, but overall, 
transition disks tend to have much lower masses and accretion rates than "full disks." 
Eight of our targets are multiples: 6 are binaries and the other 2 are triple systems. In 
four cases, the stellar companions are close enough to suspect they are responsible for the 
inferred inner holes. We do not see an increased incidence of binaries with separations in 
the ~8-20 AU range, the range where we are sensitive to companions that could carve the 
inferred inner holes, suggesting companions at these separations are not responsible for a 
large fraction of the transition disk population. However, given the small size of the current 
sample this result should not be overinterpreted. A complementary radial velocity survey to 
find the tightest companions is highly desirable to firmly establish the fraction of transition 
disks that could be accounted for by very tight binaries. 

We find that 9 of our transition disk targets have low disk mass (< 2.5 Mjup) and 
negligible accretion (< 10~^^ Moyr"^), and are thus consistent with photoevaporating (or 
photoevaporated) disks. Four of the non-accreting objects have fractional disk luminosities 
< 10^^ and could already be in the debris disk stage. The remaining 17 objects are 
accreting. Four of these accreting objects have SEDs suggesting the presence of sharp inner 
holes {aexcess values > 0), and thus are excellent candidates for harboring giant planets. 
The other 13 accreting objects have aexcess values < 0, which suggest a more or less radially 
continuous disk. These systems could be forming terrestrial planets, but their planet 
formation stage remains unconstrained by current observations. 

Understanding transition disks is key to understanding disk evolution and planet 
formation. They are systems where important disk evolution processes such as grain 
growth, photevaporation, dynamical interactions and planet formation itself are clearly 
discernablc. In the near future, detailed studies of transition disks such as sources # 11, 
21, 31, and 32, will very likely revolutionize our understanding of planet formation. In 
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particular, the Atacama Large Millimeter Array (ALMA) will have the resolution and 
sensitivity needed to image these transition disks, using both the continuum and molecular 
tracers, at ~l-3 AU resolution. Such exquisite observations will provide unprecedented 
observational constraints, much needed to distinguish among competing theories of planet 
formation. Finding promising targets for ALMA is one of the main goals of this paper, the 
first one of a series covering over 100 Spitzer-selected transition disks. 
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Table 3. Derived Properties 



# 


LOG(Acc. rate) 


Mass Disk'' 


r • b 


tun off 




LOG(L£,/L«) 




Object Type 




(Me/yr) 


{Mjup) 


(AU) 








(mag) 




1 


-9.3 


< 1.1 




2.20 


-1.16 


-1.88 


0.7 


grain growth dominated disk 


2 








8.00 


-1.25 


< -4.09 


0.2 


Be star 


3 


-8.0 


3.3 




2.20 


-0.81 


-1.68 


1.1 


grain growth dominated disk 


4 








8.00 


-2.28 


-3.92 


1.7 


AGB star 


5 


< -11.0 


< 0.8 


68 


5.80 


0.38 


-2.32 


0.7 


photoevaporating disk 


6 








8.00 


-2.01 


< -3.92 


1.3 


AGB star 


7 


-9.5 


< 1.4 




4.50 


-0.63 


-2.55 


0.4 


grain growth dominated disk 


8 








8.00 


-2.00 


< -3.70 


1.4 


AGB star 


9 


< -11.0 


< 1.5 


210 


2.20 


-1.20 


-1.50 


0.9 


photoevaporating disk 


10 








8.00 


-2.14 


< -3.72 


2.0 


AGB star 


11 


-8.8 


< 1.5 




8.00 


0.65 


< -2.72 


1.0 


giant planet forming disk 


12 


< -11.0 


< 1.1 


0.62 


8.00 


-0.21 


-3.96 


2.0 


circumbinary/debris disk 


13 


< -11.0 


< 0.6 




5.80 


-0.99 


-2.50 


0.1 


photoevaporating disk 


14 


-7.2 


37.5 




2.20 


-1.25 


-2.06 


1.2 


grain growth dominated disk 


15 


-9.6 


0.6 


69 


4.50 


-1.27 


-2.51 


2.7 


grain growth dominated disk 


16 


-9.7 


< 0.7 




2.20 


-1.25 


-2.01 


1.9 


grain growth dominated disk 


17 


-9.6 


5.4 




4.50 


-0.95 


-2.44 


1.4 


grain growth dominated disk 


18 


-9.4 


1.2 


81 


2.20 


-1.02 


-1.78 


0.8 


grain growth dominated disk 


19 


< -11.0 


< 1.3 




8.00 


0.16 


< -2.96 


1.1 


photoevaporating disk 


20 


< -11.0 


< 1.6 




8.00 


-1.99 


< -3.21 


1.6 


debris disk 


21 


-9.3 


1.3 




8.00 


0.69 


-2.30 


1.4 


giant planet-forming disk 


22 








8.00 


-1.92 


< -3.03 


2.6 


AGB star 


23 


< -11.0 


<1.7 


24 


8.00 


-2.18 


< -3.76 


1.0 


circumbinary/debris disk 


24 


< -11.0 


<2.1 


105 


5.80 


-1.07 


-2.85 


0.0 


circumbinajry/photoeva. disk 


25 


< -11.0 


< 0.7 




5.80 


-1.56 


-3.41 


0.7 


debris disk 


26 


< -9.7 


< 1.2 




4.50 


-1.00 


-2.44 


0.6 


grain growth dominated disk 


27 


-8.4 


< 0.8 


41 


4.50 


-0.62 


-2.38 


0.7 


circum. /grain growth dominated disk 


28 


-9.3 


< 2.1 




2.20 


-1.49 


-2.06 


1.3 


grain growth dominated disk 


29 


-8.8 


< 2.0 




3.60 


-1.31 


-2.24 


1.6 


grain growth dominated disk 


30 


-8.3 


< 1.1 




4.50 


-0.21 


-2.23 


2.6 


grain growth dominated disk 


31 


-7.3 


< 1.3 




8.00 


0.30 


< -2.73 


1.3 


giant planet forming disk 
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Table 3 — Continued 



# 


LOG(Acc. rate) 


Mass Disk'' 


r 

'-pro J . 


^tun — off 


O^excess 


LOG(Lb/L,) 


Aj 


Object Type 




(Mo/yr) 


(Mjup) 


(AU) 


fim 






(mag) 




32 


-9.9 


11.1 




8.00 


0.72 


-2.70 


1.3 


giant planet forming disk 


33 








8.00 


-2.21 


< -3.86 


1.6 


AGB star 


34 








8.00 


-1.83 


< -3.58 


1.4 


AGB star 



^The disk mass upper limits for targets #12 and 27 come from Cieza et al. (2008). 

''Source # 12 is a binary identified by VLBA observations (Loinard et al. 2008). Source # 24 is a triple system. The tight 
components are consistent with two equal-brightness objects with a separation of ~7 AU and a ^30 deg position angle (see 
§ |4.4| and Figure [4]|. Source # 27 is a triple system. The "primary" star in the VLT observations is itself a spectroscopic 
binary with a 35.9 d period and estimated separation of 0.27 AU (Mathieu et al. 1994). 
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Fig. 1. — Color-color diagram illustrating our key target selection criteria. Objects with 
[3.6]-[4.5] < 0.25 and [3.6] -[24] < 0.5 are consistent with bare stellar photospheres. Blue 
dots are disk-less WTTSs from Cieza et al. (2007) used to define this region of the diagram. 
Red stars are all the 297 Young Stellar Objects Candidates (YSOc) in the Cores to Disks 
catalog of Ophiuchus. Most PMS stars are either disk-less or have excesses at both 4.5 and 
24 fim. Our 26 transition disks, shown as green dots, have significant (> 5-a) excess at 24 
/im and little or no excess at 4.5 fim, as expected for disks with inner holes. The 8 black 
dots are likely Be and AGB stars contaminating our original sample. 
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-200 200 -200 200 

Velocity (km/s) 



Fig. 2. — The Ha velocity profiles of the 17 accreting objects in our sample. The dashed line 

indicates the 10% peak intensity, where AV is measured. They all have AV > 300 km/s. 
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Fig. 3. — The Ha velocity profiles of the 8 non-accreting objects in our sample where Ha was 
detected. The dashed hne indicates the 10% peak intensity, where AV is measured. Non- 
accreting objects show symmetric and narrow (AV < 230 km/s) Ha emission, consistent 
with chromospheric activity. 
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Fig. 4. — The K-band images of the seven multiple systems that have been detected with our 
VLT-AO observations. Target #24 is a triple system. The tighter components (#24-B and 
#24-C) are not fully resolved, but they presence can be inferred from the highly elongated 
image less than 1" away from the perfectly round PSF of the source ^2A-A 
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Fig. 5. — The [3.6] -[24] vs [3.6]-[4.5] colors of Ophiuchus PMS stars with disks of different 
masses (left) and accretion rates (right). Transition disks (in the left side of the figure, 
with [3.6]-[4.5] < 0.25 ) tend to have much lower disk masses and accretion rates than non- 
transition objects. The disk masses of non-transition disk sample have been taken from 
Andrews & Williams (2007). The accretion rates of non-transition disks have been taken 
from Natta et al. (2006). 
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Fig. 6. — Disk masses vs. accretion rates for accreting objects with mm detections (blue 
circles), accreting objects with mm upper limits (red arrows), and non-accreting objects with 
mm upper limits (green triangles). None of the non-accreting objects was detected at mm 
wavelength. 
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Fig. 7. — Distribution of companion separations for our transition disk sample (red dashed 
line), non-transition disks (blue dash-doted line), and disk-less stars (green dotted line). 
Spectroscopic binaries have been assigned a projected separation of 0.5 AU. The total number 
of objects (i.e., single + multiple systems) in each sample are shown in parenthesis. The 
distribution of binarv senarations for main seauence fMS") solar-tvne stars fDunuennov & 
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Fig. 8. — dexcess VS. Xturn-off for objects with different L^/L* values. The loci of "typical" 
CTTS with full disks are shown for comparison. Transition disks occupy a much larger 
parameter space 



- 58 - 





-8 




-10 




-12 


1 




a 










-8 


1 

m 






-10 








er 


-12 






< 


-14 








-8 




iafl 




O 






-10 




-12 




-14 



IRAN 5 " 

^'^^ v; ■■• ■ 
: M„ < O.QMjup \ 


TRAN 9 " 


" TRAN 12 " 
. Mb < l.lMjup \ ^ : 


IRAN 13 " 


TRAN 19 " 
. Mo < 1.3Mjup \^ ^. 


TRAN 20 ; 
. Mb < 2.2Mjup \ X 


IRAN 23 " 

: M„ < 2.4Mjup \ ;i; 


TRAN 24 " 
: Mb < 2.1Mji,p \ X 


TRAN 25 : 
. Mb < O.VMjup \^ 



012301230123 

log A (yLtm) 



Fig. 9. — The SEDs of our non-accreting transition disk targets. All of them have low 
disk masses (< 2.5 Mjup), and are consistent with photoevaporation. The filled circles are 
detections while the arrows represent 3-0" limits. The open squares correspond to the observed 
optical and near-IR fluxes before being corrected for extinction using the Aj values listed in 
Table 3 and the extinction curve provided by the Asiago database of photometric systems 
(Fiorucci & Munari 2003). For each object, the average of the two R-band magnitudes 
(from the USNO-Bl catalog) listed in Table 1 has been used. The solid blue line represent 
the stellar photosphere normalized to the extinction-corrected J-band. The dotted lines 
correspond to the median mid-IR SED of K5-M2 CTTSs calculated by Furlan et al. (2006). 
The dashed lines are the quartiles. 
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Fig. 10. — The a.excess (the slope of the IR excess) vs. Xtum-off (the wavelength at which 
the excess becomes significant) for our sample of transition disks. The masses and accretion 
rates are indicated by different symbols. The 4 accreting (targets #11, 2, 31, and 32) and 
the 2 non-accreting (targets #5 and 19) objects with Uexcess ^ are suggestive of sharp 
inner holes. 
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Fig. 11. — The SEDs of all our accreting transition disk targets. The symbols are the same 
as in Fig. 9. Disk masses range from < 0.7 to ~40 Mjup. For targets # 14, 17, and 32, 850 
um fliixes. from Andrews & Williams f2007') and Nutter et al. flQQB). are shown in addition 



